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Optical Monitoring of PKS 1510-089: A Binary Black Hole 
System? 
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ABSTRACT 

Three deep flux minima were observed with nearly the same time-scales and intervals for the 
blazar PKS 1510—089 in the past few years. A binary black hole system was proposed to be 
at the nucleus of this object, and a new minimum was predicted to occur in 2002 March. We 
monitored this source with a 60/90 cm Schmidt telescope from 2002 February to April. In 
combination with the data obtained bv lXieetal] ll2004) in the same period, we presented for 
the 2002 minimum a nearly symmetric light curve, which would be required by an eclipsing 
model of a binary black hole system. We also constrained the time-scale of the minimum to 
be 35 min, which is more consistent with the time-scales (~ 42 min) of the three previous 
minima than the 89 min time-scale given by the same authors. The wiggling miniarcsecond 
radio jet observed in this object is taken as a further evidence for the binary black hole system. 
The 'coupling' of the periodicity in light curve and the helicity in radio jet is discussed in the 
framework of a binary black hole system. 
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1 INTRODUCTION 

Blazars, as a subset of active galactic nuclei (AGN), are charac- 
terized by rapid and strong variability in multi-wavebands. This 
behaviour gives us much important information on their central 
physics. Variability studies of blazars have been essential in un- 
derstanding their inner structures, radiation mechanism, and other 
physical processes. 

PKS 1510—089 has been the target of m any monitoring pro- 
grams (e.g., Villata et al. 1997; Xie et al. 2002). It is a flat-spectrum 
radio quasar at redshift 0.361 and has a pars ec-scale jet wi thin 
just 3° to our line of sight according to iHoman et alj 120021) . A 
pronounced UV-excess, a very flat X-ray spectrum, and a steep 
gamma-ray spectrum are found in this object. Its light curve is quite 
different from those of other blazars. Besides fast and small am- 
plitude variations and irregular outbursts, some brightness minima 
were also observed (see lXie et al^OO^ and ref erences there in) . 
In particular, three deep minima were reported by Dai et al. 1 2001 ), 
IXieetalJl200ll.l2002t) with variations of 0.65 mag/41 min on 1999 
June 14, 2.00 mag/41 min on 2000 May 29, and 0.85 mag/44 min 
on 2001 April 16, respectively. These minima have nearly the same 
time-scales (~ 42 min) a nd intervals (336 ± 14 d), and based 
on which. Ixie et alj {2002) argued that the minimum was a peri- 
odic phenomenon and proposed that there was a binary black hole 
(BBH) system in the centre of PKS 1510—089. The minimum oc- 
curs when the primary black hole (PBH) is eclipsed by the sec- 
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ondary black hole (SBH). They also predicted that a new minimum 
would o ccur in 20 02 March, and they confirmed it with new obser- 
vations IXie et all2004l) . 

However, t he time-scale of the new minimum reported by 
IXie et al] 120041) is 89 min, which is more than two times of the 
time-scales of the three previous minima. In fact, the authors had no 
observation for about 55 min before the minimum point (see figs. 8 
and 9 in lXie etaill2004h . making their recorded fading phase (73 
min) much longer than the brightening phase (16 min). This highly 
asymmetric light curve would be rejected by an eclipsing model of 
a binary system. 

We also monitored this obj ect from 2002 February to April. 
In combination with the data in lXie et al. (2004), we provided new 
constraints on the optical variations of PKS 1510—089 in this pe- 
riod. Sect. 2 describes our observation and data reduction proce- 
dures. The results are presented in Sect. 3. Sect. 4 discusses the 
evidence for BBH from periodic light curves and wiggling jets and 
Sect. 5 gives a summary. 



2 OBSERVATION AND DATA REDUCTION 

Our optical observation was performed on a 60/90 cm Schmidt tele- 
scope located at the Xinglong Station of the National Astronomical 
Observatories of China (NAOC). A Ford Aerospace 2048 x 2048 
CCD camera is mounted at its main focus. The CCD has a pixel 
size of 15 fi and a field of view of 58' x 58', resulting in a 
resolution of l'.'7 pixel - . The telescope is equipped with a 15 
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Figure 1. Finding chart of PKS 15 10— 089 and four comparison stars used 
by us. The image was taken with our Schmidt telescope in the BATC i band 
on JD 2 452 344 with size 8' X 8' (cut from a 58' X 58' image). 
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Figure 2. Light curve of PKS 15 10-089 in the R band from 2002 Feb 23 
to Apr 10. 



color intermediate-band photometric system, which covers a wave- 
length range from 3000 to 10 000 A. The telescope and photometric 
system are mainly used to carry out the Beijing-Arizona-Taiwan- 
Connecticut (BATC) survey and have shown their efficie ncy in de- 
tectin g fast variability in blazars (e.g., Peng et al. 2003; Wu et alJ 
2005). 

Our monitoring program covered the period from 2002 Febru- 
ary 23 to April 10. Three or four photometric measurements were 
made in each night. However, as a result of the weather condition, 
only 14 night's data are useful. We used the most sensitive i fil- 
ter of the BATC photometric system, which has a central wave- 
length of 6711 A and a passband of 497 A, and which is close to 
the Cousins' R band. An exposure time of 300 ~ 600 s, depending 
on the weather and seeing conditions, was able to produce a CCD 
image with a good signal-to-noise ratio. The observational parame- 
ters are summarized in TableQ Fig.Qillustrates the finding chart of 
PKS 1510—089 and the four comparison stars used by us. The four 
comparison stars were also used in other monitoring program (e.g., 
Villat a et alJll997tlx7e et alJl2002h and have the R magnitudes of 
13.95 ± 0.03, 14.22 ± 0.03, 14.35 ± 0.05, and 14.61 ± 0.02, re- 
spectively. 

Two automatic procedures have been developed for the data 
reduction of the BATC images. 'Pipeline F includes bias subtrac- 
tion and flat-fielding of the CCD images, and 'Pipeline IF mea- 
sures the instrumental magnitudes of point sources in the BATC 
images. The latter is base d on Stetson's standard procedure of 
DAOPHOT <Stetsonlll987T) . The extinction coefficients and zero 
points o f the inst rumental magnitudes were obtained by observing 
the four lOke & Gunnl ll983l) standard stars and were used to cali- 
brate the instrumental magnitudes into the BATC AB magnitudes. 
For details of the data reduction procedures, see IZhouetaljl2003l) . 

The light curv es of PKS 15 10— 089 were given in the R band 
in IXie et"all 120041) . For comparison, we calibrated our i magni- 
tudes into the R magnitudes, which was done by using the perfect 
linear relationship between the BATC i magnitude and the Cousins' 
R magnitude found bv lZhouetaTl <2003l) . i.e., R = i + 0.1. Out- 
light curves in the R band were then obtained. 



3 LIGHT CURVES AND TIME-SCALE OF THE 
MINIMUM 

The observational results are given in TableQwith columns being 
the observational date and time (UT), Julian Date, exposure time, 
and the R magnitude and its error. The overall light curve in the R 
band is plotted in Fig.|2| in which data in IXie etalJ f2004) are also 
plotted. The two data sets are basically consistent with each other, 
except that our amplitudes of variation within each night are much 
smaller than those of Xie et al.'s, which is possibly due to our much 
shorter monitoring durations within each night. 

The predicted 2002 minimum occurred on March 8 
(JD 2452 342) according to IXie et alJ <2004 . The light curve in 
that night is given in Fig. [3] which combines our and Xie et al.'s 
data. The upper panel illustrates the light curve of the target quasar 
while the lower panel shows the differential magnitude (average 
set to 0.0) between the 2nd and 4th comparison stars. From Fig. [3] 
we can see that our observations just filled the large gap of Xie et 
al.'s observations. After a short time of fading phase beginning at 
JD 2 452 342.329, the source underwent a brief brightening phase, 
as observed by us but missed by Xie et al. Then its brightness in 
the R band dropped from 16.62 ±0.12 to 17.92 ±0.15 mag within 
19 min, and rose back to 16.57 ± 0.15 mag in 16 min. Therefore, 
the actual ti me-scale of the m inimum is 35 min, rather than 89 min 
reported bv lXie et alJ 120041) . The light curve of the minimum be- 
comes basically symmetric, in contrary to Xie et al.'s result. 

As mentioned in Sect. 1. iDai et alJ feOOlh . Ixle et alJ <200lL 

2002) reported th ree minima with time-scales of 41, 41, and 44 
min, respectively. IXie et al] 120021) argued that the minimum is a 
periodic phenomenon and there is a BBH system in the centre of 
PKS 1510—089. They used the measured time-scales to estimate 
the black hole masses. In this sense, the new 35 min time-scale of 
the 2002 minimum, as constrained by a combination of our and 
IXie et alJ Jiool data, is more consistent with the time-scales of 
the th ree previou s minima than the 89 min time-scale measured by 
IXie et al] 12004). Moreover, the basically symmetric light curve is 
more reasonable for an eclipsing model. 
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Figure 3. Light curve of PKS 1510-089 in the R band on 2002 Mar 8. 
The open circles and dotted line are data in lXie et afl 1200 4) while the rilled 
circles and dashed line are our observations. The lower panel gives the dif- 
ferential magnitudes between the 2nd and 4th comparison stars. 



4 DISCUSSION 

4.1 BBH Model from Periodic Light Curves 

Periodicities in light cu rves are o ften attributed to the presence of 
BBHs (see a review bv lKomossiJl2003l and references therein). A 
prominent example is the BL Lac object OJ 287, which shows in 
the optical bands quite a strict period of 1 1.86 vr |Siilanpaa et alJ 
1 1 98 8l 1 1 99d I Valtaoi a et all2000tlPursimo et all200d) . Most exphT 
nations of the periodicity are based on the as sumption that there is 
a BBH system in the nucleus of OJ 287 (e.g. JSillanp aa et alJl988t 
IVillata et alll998tlValtaoia et all200d:lLiu & Wuj2002l) . The peri- 
odicity is linked to the orbital motion of the BBH system. 

For PKS 1510-089, the situation is much different. Its light 
curves show periodic minima rather than periodic outbursts. The 
best model to explain the mini ma in such short time-scales is the 
eclipsing model, as described bv lXie et alJ<2002h . Other possibili- 
ties can hardly result in minima with such short time-scales. In fact, 
the four minima with nearly the same time-scales and intervals and 
the basically symmetric light curves would argue strongly for an 
eclipsing model and a BBH system. 



4.2 Evidence for BBH from Radio Observations 

Evidence for BBH in PKS 1510—089 comes also from radio ob- 
servations. PKS 1510—089 is a radio-loud source with high po- 
larizations. Earl y radio observations de tected an arcsecond jet in 
this object lO'Dea. Barvainis , <^Oialli^P98^) and classified it as a 
source wi th a well-aligned jet jTinga^Murphv. & Edwardsl Il998t 
ICaol20 00). More recently, the more sensitive VLBA and VLBI ob- 
servations detected a mini arcsecond (mas) jet in the o pposite di- 
rection of its arcsecond je t iFev & Charloil 9*971 [Kellerman n et alJ 
Il998t Ijorstadetal]l200lt iHoman et alj|200ll) . and then the ob- 
ject shows perhaps the most highly misaligned jet ever observed 
IHoman et all2002h . 

In order to connect th e highly misaligned mas and arcsec- 
ond iets. lHoman et alJ J2002l) argued that the jet in PKS 15 10-089 
bends at a deprojected distance of ~ 30 kpc from the nucleus, near 
the probable boundary of the host galaxy. The bend occurs either 




SBH Orbit 

Figure 4. Schema of the configuration of the BBH system in 
PKS 1510-089. The Earth is in the orbital plane of the SBH. When the 
SBH moves along the orbit, the inner part of the primary accretion disc will 
vibrate between positions Dl and D2, and the jet will oscillate between Jl 
and J2, resulting in a wiggling pattern. (The jet is not in the orbital plane and 
actually follows a highly eccentric precess when the SBH moves along its 
orbit. See text for details.) The angles are enlarged for clarity. The minimum 
occurs when the SBH moves to position B3. 



by ram pressure from winds in the intracluster medium or by the 
density gradient in the transition to the intergalactic medium. 

In fact, the jet shows bends on all scales on the 1.7 and 5 
GHz VLBA maps, from mas to arcsecond. The smaller the scale , 
the more s ignificant the bend (see fig. 2 in IHoman et al]l2002h . 
IHoman et alJ l2002h also noted that the jet ridge-line 'wiggles in 
the plane of the sky after the first 10 mas'. Similarly, the 15 GHz 
image of Keller mann et alJ ll998h reveals a short jet extending to 
the north and then bending to the northwest within 2 mas of the 
core. Furthermore, multi-epoch VLBA observations have detected 
several components moving away from t he nucleus wit h different 
speeds and position angles (see fig. 26 in ljorstad et all200iri . This 
pattern of ejection naturally leads to a wiggling or precessing inner 
jet. 

The wiggling or precessing jet can be best explained in terms 
of a BBH system. BBHs were first proposed to be in AGNs 
by Bcgclman. Blandford. & Ree| Jl980h and have since been fre- 
quentl y used to ex plain the observ ed curved or mis aligned jets (e.g., 
iRoos. Kaastra. & Hummelll993t Icoiiwav & Wrobe1ll995t) . BBH 
systems are formed mainly via galaxy mergers. In most cases, the 
SBH is in an orbit non-coplanar with the accretion disc of the PBH, 
inducing torques in the inner parts of the disc and resulting in pre- 
cession of the disc and its jet. 

Therefore, the VLBA and VLBI observations of the distorted 
jets in PKS 15 10— 089 provide another evidence for the presence of 
the BBH in its nucleus. However, the situation in PKS 1510—089 
is much different. Its light curves show periodic minima rather than 
periodic outbursts. Correspondingly, it is also different in the in- 
terpretation of the periodic minima, which is linked to the peri- 
odic e clipse of the PBH by the naked SBH according to lXie etal] 
(2002). In their scenario, the orbit of the SBH passes our line of 
sight to the PBH. This configuration represents an extreme case of 
BBH systems, for which a schema is given in Fig. |4] The circle 
represents the orbit of the SBH around the PBH, and the Earth is in 
the orbital plane. The thick lines denote two extreme positions of 
the inner part of the primary accretion disc. When the SBH moves 
to position Bl, the inner part of the primary accretion disc tends 
to move to position Dl under the gravity of the SBH, and its jet is 
directed to Jl. When the SBH moves to position B2, the inner part 
of the primary disc tends to move to position D2 under the gravity 
of the SBH, and its jet is directed to J2. Therefore, the SBH moves 
along the circle, makes the inner part of the primary disc vibrat- 
ing between Dl and D2 and the jet oscillating between Jl and J2. 
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Table 1. Observational log and results 



Observation Date 
(UT) 


Observation Time 
(UT) 


Julian Date 


Exposure Time 
(s) 


R 
(mag) 


FtcTT 

(mag) 


700? 09 91 

ZUUZ Ui. Z.J 


90-96-91 


94.5919Q 159 

ZtJZJZ7.JJZ 


300 


16 83 


04 


700? n? 91 


90-11-58 
ZU.jj.JO 


94.5919Q 157 
ZtJZjZy.JJ / 




10.01 




9009 09 91 

ZUUZ UZ ZJ 


20"41"36 


94.5919Q 169 

ZtJiJZ7.JvZ 


300 


16.83 


0.04 


9009 01 OS 


90-4.7-4.0 
ZU.t / .tU 


94.5911Q 166 

ZUZJJ7.JUU 


ftftfl 

□UU 


16 80 


OS 


9009 01 OS 


91 -00-90 

Z 1 .UU.ZU 


94.5911Q 175 

ZtJZJJ7.J / J 


600 


16 77 


0.05 


9009 01 06 
ZUUZ UJ UO 


90-14.-5Q 

ZU.Jt.J7 


94.5914.0 158 
Zt J ZJtU. JJ 




16 71 


OS 


9009 01 06 

ZAA./Z. UJ UU 


90-4.9-98 
zu.tz.zo 


94.5914.0 161 

ZtJZJtU.JUJ 


300 


16 75 


0.05 


9009 01 0ft 
ZUUZ UJ UO 


90-50- 1 7 
ZU.jU. 1 / 


94.5914.0 168 
ZtJZJtU.jUO 




16 70 


OS 


9009 01 07 

ZUUZ U J u / 


90-98- 1 6 

_w,_0, 1 U 


94.5914,1 151 

ZtJZJtl . Jj J 


300 


16 69 


05 


9009 01 07 
ZUUZ UJ U / 


90-15-4.8 
ZU.JJ .'-to 


94.5914.1 158 
ZtJZj'tJ .JJO 






05 


9009 01 07 

ZUUZ U J u / 


20"43"40 


94.5914.1 164. 

ZtJLJt 1 .JUt 


300 


16.66 


05 


9009 01 07 
ZUUZ Uj U / 


90-5 1 -OR 

ZU. J 1 .UO 


94.5914.1 16Q 

ZUZj^ 1 .JU7 




16 78 


06 
u.uo 


9009 01 OS 

Z.UV7Z. UJ UO 


90-11-04. 

ZU. J J .ut 


94.5914,9 15ft 

ZtJ ZJtZ .JJU 


300 


16.47 


0.12 


9009 01 08 
ZUUZ Uj UO 


90-4.0-4.7 
zu.tu.t / 


94.5914.9 169 
ZtJZjH-Z.jUZ 




16 60 


OQ 


9009 01 08 

Z.UV7Z. UJ UO 


20"48"20 


94.5914.9 167 

ZtJZJtZ.JU / 


300 


16.62 


0.12 


9009 01 0Q 
z.uuz. uj \jy 


20-45-13 


94.5914.1 165 

ZtJZJtJ.JUJ 


300 


16 75 


0.03 


9009 01 0Q 

ZUUZ UJ U7 


21 "01 "52 


9459141 176 

ZtJ ZJtJ . J / u 


300 


16 71 


0.03 


9009 01 0Q 

Z.UVZ UJ U7 


21"09"24 


94.5914,1 189 

ZtJZJtJ . JOZ 


300 


16.65 


03 


2002 03 1 

_ i iVJ— U J A U 


90-1 ,-07 


2452344 358 

itJ iJTT.JJ O 


300 


16.69 


0.04 


9009 01 1 

Z.UV7Z. UJ 1 u 


20"42"37 


94.5914.4. 161 

ZtJZJtt.JUJ 


300 


16 68 


04 


9009 01 1 

LUUZ U J 1 u 


90-50-0Q 

ZU.JU.U7 


94.5914,4. 168 

ZtJZJtt.JUO 


300 


16.67 


04 


9009 01 1 

ZUUZ UJ 1 u 


90-57-50 

ZU.J / . JU 


94.5914,4. 171 

ZtJZJtt.J / J 


300 


16 69 


0.04 


9009 01 1 9 

Z.UUZ. U J 1 z 


20-07-37 


94.5914.6 HQ 

ZtJLJtU.JJ" 


600 


16.66 


0.05 


9009 01 1 9 

Z.UUZ. U J 1 z 


90-90-1 7 
zu.zu. 1 / 


94.5914.6 14.7 

ZtJiJtU.Jt / 


600 


16.65 


0.04 


9009 01 14. 

Z.UUZ UJ It 


1 7.JU.UU 


94.5914,8 111 

ZtJZJtO.JJ 1 


600 


16.61 


0.02 


9009 01 14. 

Z.UUZ UJ It 


90-08-4.5 

ZU.UO.tJ 


94.5914,8 HQ 

ZtJZJtO.JJ7 


600 


16.67 


0.02 


9009 01 1 4. 

Z.UUZ. UJ it 


20"21"29 


94.5914.8 14.8 

ZtJZ.JtO.JtO 


600 


16.62 


0.02 


9009 01 91 

Z.UUZ. UJ ZJ 


19"30"26 


94.59157 111 


600 


16 69 


03 


9009 01 91 

Z.UUZ. UJ ZJ 


19-43-47 


94.59157 199 

ZtJZJJ / . JZZ 


600 


16.72 


03 


2002 03 23 


19:56:18 


2452357.331 


600 


16.66 


0.03 


2002 03 24 


20:32:04 


2452358.356 


300 


16.70 


0.05 


2002 03 24 


20:40:39 


2452358.362 


300 


16.58 


0.05 


2002 03 24 


20:48:24 


2452358.367 


300 


16.59 


0.05 


2002 03 26 


20:31:41 


2452360.355 


600 


16.74 


0.06 


2002 03 26 


20:44:20 


2452360.364 


600 


16.76 


0.09 


2002 03 26 


20:56:59 


2452360.373 


600 


16.68 


0.07 


2002 04 09 


19:38:24 


2452374.318 


600 


16.78 


0.02 


2002 04 09 


19:50:52 


2452374.327 


600 


16.77 


0.02 


2002 04 09 


20:03:30 


2452374.336 


600 


16.83 


0.03 


2002 04 10 


19:46:14 


2452375.324 


600 


16.70 


0.05 


2002 04 10 


19:58:51 


2452375.333 


600 


16.81 


0.03 


2002 04 10 


20:11:27 


2452375.341 


600 


16.86 


0.03 



Consequently, an apparently wiggling jet is formed. When the SBH 
moves to B3, the point on our line of sight to the PBH, the PBH is 
eclipsed and a flux minimum occurs. The minimum is a periodic 
phenomenon, and the periodicity is caused by the orbital motion of 
the BBH system. 

For simplicity, Fig. |4] gives only a 2-D diagram. In the 3-D 
reality, the jet will have a small angle to the orbital plane of the 
SBH, and the primary accretion disc will be slightly tilted relative 
to the normal of the plane. When the SBH moves along its orbit, 
the jet and the primary accretion disc will follow a highly eccentric 
precess, with Jl and J2 being two extreme directions of the jet. The 
angles in Fig. [4] are enlarged for clarity. The jet actually stays in 
roughly the same direction from the point of view of the observer 
with only a small wiggling in the plane of the sky. During the orbital 
motion of the BBH system, the SBH could eclipse the PBH but it 
does not pass through the jet, and the jet does not cross the line of 
sight when precessing. 



4.3 The 'Coupling' of Periodicity and Helicity in Blazars 

It is interesting to notice that (quasi-)helical jet and pe- 
riodic light curve usually co-exist in blazars. OJ 287 it- 
self s hows evidence for a helical jet in its 8.4 GHz VLBI 
map 1 Vicente. Chariot. & Soli Il996i), Other examples present- 
ing bo t h helic a l jets a nd perio dic light curves inclu de 3C 120 
jWebbl Il99d: iGomez. Marscher. & Alberdil Il999bh. 3C 34 5 
iwebb. Smith. & Leacocklll988l: IZensus. Cohen & Unwinlll995h . 
AO 0235+16 llorstad et all I200TI: Raiteri et all 1200 lll . BL Lac 
jTatevama et all Il998t iFan et all 1 19981) . Mrk 501 fai et all 
I l995t lHavashida et all Il998l) . PKS 0735+17 iFan et all 119971 
Gome z et all Il999al) . and these objects are most proposed to 
harbor a BBH system in their nuclei {V ill ata & Raiteril T 



ICaproni & Abraham! f2004albl lOstorero. Villata. & Raiteril I 
In other words, the periodicity in light curve and helicity in jet usu- 
ally 'couple' in blazars and are interpreted within the scenario of a 
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BBH system. In fact, the 'coupling' of the periodicity and helicity 
can be easily explained in that the periodicity in helical jets will 
result in periodicity in relativistic boosting in jets and hence the 
periodicity in brightness. 

Similarly, the periodic minima and wiggling jet observed in 
PKS 1510—089 can both be taken as evidence for a BBH system. 
The SBH orbits around the PBH and exerts tidal torques on its ac- 
cretion disc, resulting in disc vibration and jet wiggling. The eclipse 
of the PBH by the SBH leads to the periodic flux minima. 



5 SUMMARY 

The blazar PKS 1510—089 was monitored from 2002 F ebruary to 
April. Based on a combination of our data and those of IXie et al] 
(2004), we constrained a new time-scale of 35 min and a nearly 
symmetric light curve for the 2002 minimum reported by Xie et al. 
The new res ults are more con s istent wit h the three previous minima 
observed bv lDai et alJ <200lMXie eta?] l200lll2002l) and are more 
reasonabl e for a n eclipsing model of a BBH system suggested by 
IXieetalj f2002). Radio observations have revealed wiggling jets in 
PKS 1510—089 and provided another evidence for the BBH sys- 
tem. The 'coupling' of the periodicity in light curve and helicity in 
radio jet is discussed within the framework of BBH systems. 
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